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INTRODUCTION 

High  costs  and  injuries  to  test  subjects  during  emergency  evacuation 
demonstrations  have  led  to  the  development  of  computer  simulation  models  of 
emergency  evacuations  with  a long-range  goal  of  eventually  replacing  the 
certification  demonstrations.  Besides  the  obvious  advantages  of  low  cost 
and  no  human  injury,  a computer  model  can  be  used  to  simulate  a wide  range 
of  emergencies  rather  than  those  done  under  the  restrictive  criteria  estab- 
lished for  a certification  demonstration. 

An  emergency  evacuation  computer  model  described  in  Reference  1 has 
been  developed  by  FAA  Aeronautical  Center  personnel.  This  computer  program 
is  based  on  the  assumption  that  the  door  is  a bottleneck  and  passenger  move- 
ment in  the  interior  of  the  aircraft  is  not  considered  except  as  a queue  at 
the  exit  door.  The  computer  model  in  Reference  1 is  written  in  GPSS . The 
present  computer  model  is  written  in  FORTRAN  and  traces  a passenger  from  his 
seat  to  the  exit  he  selects  during  evacuation.  Separate  programs  have  been 
developed  for  wide  and  narrow  body  aircraft. 


MODEL  DESCRIPTION 

AIRCRAFT  INTERIOR  REPRESENTATION  - The  interior  of  the  aircraft  is  treated  as 
a matrix  of  I,  J values  that  identify  passenger  seating,  aisle,  and  exit 
locations.  A general  representation  of  wide  and  narrow  body  aircraft  is 
shown  in  figures  1 and  2.  The  I denotes  row  number  and  the  J denotes  column 
number.  Thus,  1*1,  J « 1 represents  the  seat  of  the  upper  left  hand  pas- 
senger in  figures  1 and  2.  Figures  1 and  2 show  input  parameters,  which 
will  be  discussed  in  detail  in  the  program  input  section,  that  allow  the 
program  user  to  specify  the  seating  arrangement  and  passenger  assigned  to  a 
particular  exit.  Further  flexibility  is  added  to  the  model  by  allowing  the 
user  to  identify  empty  seats  in  the  aircraft.  In  the  narrow  body  program 
the  user  can  identify  seats  in  exit  rows  that  are  common  in  overwing  exits. 

EXIT  PATH  TIMS  COMPONENTS  - Referring  to  figures  1 and  2,  in  evacuating  a 
passenger  must  get  out  of  his  seat  and  move  through  boxes  where  seats  are 
present,  move  through  aisle  boxes,  and  then  move  through  boxes  which  may  or 
may  not  contain  seats,  to  reach  the  exit  door.  Let  T^  be  the  time  required 
for  a passenger  to  move  through  a box  without  seats  and  T2  the  time  required 
for  a passenger  to  move  through  a box  with  seats. 

As  a passenger  moves  toward  the  exit  door  he  may  encounter  a delay  due 
to  a line  at  the  door.  This  delay  will  be  dependent  on  the  movement  of  the 
passengers  in  front  of  him  and  the  door  opening  time.  Let  T3  be  the  time 
required  for  door  opening  and  exit  equipment  deployment. 

When  a door  exit  is  free,  the  passenger  goes  through  the  door  and  down 
a slide  or  off  a wing  to  reach  the  ground.  Let  T4  be  the  time  required  for 
a passenger  to  go  through  the  exit  door  and  T^  the  time  required  for  a 
passenger  to  reach  the  ground. 
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The  total  evacuation  time  of  a passenger  can  be  represented  as  a func- 
tion of  T.,  T_,  T3,  T4,  and  T$.  First  the  individual  time  components  must  be 
determined  ana  then  be  combined  in  a suitable  manner  to  determine  passenger 
evacuation  time.  The  computer  programs  have  been  written  to  accept  time  com- 
ponent input  for  each  open  exit.  All  times  except  T3  are  passenger  related. 


DETERMINATION  OF  EXIT  PATH  TIMES  - The  same  mathematical  model  is  used  for 
all  five  time  components.  Each  passenger  is  given  a unique  Tj_,  T2,  T4  and 
T- . Each  open  exit  is  given  a unique  T3.  For  each  time  component , evacua- 
tion test  data  suitable  to  that  segment  of  the  evacuation  path  are  required. 

The  mean,  standard  deviation,  minimum,  and  maximum  component  times  are  required 
for  a group  of  passengers  or  exit  doors  as  program  input.  A gamma  function 
is  fitted  to  this  data  by  the  Method  of  Moments  (2) . From  the  gamma  function 
a table  of  probability  versus  time  is  generated.  A random  probability  be- 
tween 0 and  1 is  then  determined.  The  random  probability  is  then  used  to 
generate  a time  from  the  table  of  probability  versus  time. 


Factors  influencing  the  time  components  are 


1.  T^  - aisle  width,  aircraft  attitude,  aisle  blockage,  visibility 
passenger  physical  condition 


seat  configuration,  seat  belt  removal  time,  aircraft  attitude 
passenger  physical  condition 


type  of  path  followed  to  ground,  passenger  physical  condition 


This  list  is  not  necessarily  complete.  Passenger  mental  attitude  is  certainly 
an  important  factor.  A systematic  small  scale  test  program  will  probably  be 
required  to  generate  suitable  input  data  for  these  functions.  The  factors 
cited  above  must  be  considered  in  obtaining  ranges  of  data. 


Because  the  data  used  to  generate  the  time  components  are  statistical,  provi' 
sion  is  made  in  the  program  to  allow  the  user  to  generate  a desired  number 
of  emergency  evacuations.  Average  data  obtained  from  a series  of  runs  is 
probably  more  representative  of  evacuation  tests.  Generally,  similar  evacua' 
tion  tests  will  produce  different  results. 


CALCULATION  OF  EVACUATION  TIME 
ger  to  a given  open  exit  door, 
time  to  the  door  (T_)  is: 


where  « number  of  boxes  without  seats  passenger  must  move  through 
to  reach  exit 

N2  - number  of  boxes  with  seats  passenger  must  move  through  to 
reach  exit 


Upon  reaching  the  door  the  passenger  is  either  delayed  by  the  time  it 
takes  the  door  to  open  or  he  goes  through  the  door.  The  time  spent  inside 
the  aircraft  (TA)  is  taken  as  the  greater  of  T3  or  TD.  The  total  evacuation 
time  (T)  is: 


T * TA  + T4  + T5  2 

If  there  is  no  delay,  the  second  passenger  to  exit  reaches  the  door  in 
the  time  determined  from  equation  1.  Note  that  N^,  T^,  N2,  and  T2  are  dif- 
ferent for  the  second  passenger.  If  the  second  passenger  is  delayed,  his 
time  to  the  door  is  given  by: 


where  and  T^  are  values  for  the  first  passenger  preceding  him.  The  T^ 
for  the  second  passenger  is  taken  as  the  greater  TD  computed  by  equations 
1 and  3.  The  total  evacuation  time  for  the  second  passenger  is  given  by 
equation  2 using  T4  and  T-  appropriate  to  the  second  passenger.  For  each 
passenger,  TA  is  used  to  keep  a record  of  the  delays  encountered  by  that 
passenger  or  is  the  time  determined  by  equation  1 depending  on  which  is 
greater.  Equation  2 is  used  to  compute  evacuation  time  for  each  passenger. 
Thus,  the  door  may  not  necessarily  be  a bottleneck  throughout  the  evacuation 
process  in  this  computer  model. 


INITIAL  PASSENGER  EXIT  SELECTION  - Two  options  are  available  in  the  computer 
programs.  In  the  first  option,  a passenger  evacuates  through  the  nearest 
open  exit.  The  nearest  exit  is  determined  by  counting  the  total  number  of 
boxes  the  passenger  must  move  through  to  reach  each  open  exit  and  selecting 
the  exit  with  the  smallest  total  of  boxes.  In  the  second  option,  passengers 
assigned  to  a given  exit  must  evacuate  through  it  if  it  is  open.  Passengers 
assigned  to  blocked  exits  evacuate  through  the  nearest  open  exit  in  both 
options.  The  second  option  has  been  included  in  the  computer  model  because 
passengers  may  not  necessarily  evacuate  through  the  exit  nearest  them.  This 
gives  the  program  user  some  flexibility  in  selecting  an  exit  for  a passenger 
if  this  has  been  observed  in  evacuation  tests.  In  both  options,  the  passen- 
ger exit  selection  is  an  initial  one  and  passengers  at  the  ends  of  lines  of 
exits  with  long  evacuation  times  are  redistributed  to  exits  with  shorter 
evacuation  times.  The  redistribution  is  based  on  exit  flow  rates  determined 
from  the  initial  passenger  exit  selection. 

UPPER  DECK  WIDE  BODY  MODEL  - Since  wide  body  aircraft  may  have  upper  or 
lower  decks  in  addition  to  the  main  deck,  provision  has  been  added  in  the 
wide  body  program  to  handle  this  situation.  The  upper  deck  model  is  shown 
in  figure  3.  It  consists  of  a passenger  seating  configuration  with  one 
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aisle  only,  a staircase  to  the  main  deck,  and  one  exit  to  the  outside.  By 
using  I,  J notation,  the  location  of  the  staircase  and  exit  may  be  varied. 

It  should  be  noted  that  the  same  I,  J values  can  be  used  for  upper  and  main 
deck  passengers  since  the  program  handles  the  decks  separately.  The  program 
user  can  assign  any  combination  of  passengers  to  the  exit  to  the  outside  and 
staircase.  For  passengers  assigned  to  the  staircase,  the  user  must  also 
select  their  escape  exit  on  the  main  deck.  The  escape  path  and  time  segments 
for  passengers  exiting  directly  are  handled  the  same  as  lower  deck  exits 
using  time  segment  data  appropriate  to  the  upper  deck  exit.  However,  for 
passengers  using  the  staircase,  a new  model  is  needed. 

The  time  it  takes  a passenger  to  move  to  the  stair  entrance  is  similar 
to  equation  1 using  Tj^  and  T2  appropriate  to  the  upper  deck  and  the  boxes  a 
passenger  oust  move  through  to  reach  the  stair  entrance.  A new  time  function 
Tg  must  be  defined  for  passenger  movement  on  the  staircase.  This  function  is 
the  same  mathematically  as  the  five  time  functions  previously  discussed. 

Upon  reaching  the  main  deck  a passenger  moves  through  a number  of  boxes  to 
reach  the  exit.  This  time  segment  is  T2  appropriate  to  the  main  deck  exit 
times  the  number  of  boxes  the  passenger  moves  through.  The  passenger  time 
to  the  exit  is  thus  composed  of  his  time  to  the  staircase  entrance,  time  on 
the  staircase,  and  time  from  staircase  to  exit  door. 

PASSENGER  EXIT  REDISTRIBUTION  - The  aircraft  exits  must  be  numbered  as  shown 
in  figure  4 in  order  for  the  exit  redistribution  logic  to  work  properly. 

After  the  initial  evacuation,  the  exit  with  the  greatest  evacuation  time  is 
determined.  Passengers  at  the  end  of  this  exit  line  are  redistributed  to 
another  exit  in  the  following  manner: 

1.  The  exit  with  the  shortest  evacuation  time  in  the  vicinity  of  the 
exit  with  the  longest  time  is  determined. 

2.  Vicinity  of  an  exit  is  defined  as  any  exit  one  exit  away  on  either 
side  of  the  aircraft.  For  example,  referring  to  figure  4,  exits 
in  the  vicinity  of  exit  1 are  exits  2,  3,  and  4.  Exits  in  the 
vicinity  of  exit  3 are  exits  1,  2,  4,  5 and  6. 

3.  The  flow  rates  for  the  longest  exit  time  and  shortest  exit  time  in 
the  vicinity  are  computed  by: 


\ ■ <tl  - 

h - n,  - %>/ss 

where  R^  » flow  rate  for  exit  with  longest  time 
R_  * flow  rate  for  exit  with  shortest  time 


longest  evacuation  time 


T.  - exit  preparation  time  for  exit  with  longest  time 


T “ exit  preparation  time  for  exit  with  shortest  time 
3S 

Tg  * shortest  evacuation  time 

N^  3 number  of  passengers  evacuating  through  exit  with 
longest  time 

Ng  = number  of  passengers  evacuating  through  exit  with 
shortest  time 

4.  The  following  values  are  determined  by: 

tav  ■ (tl  * V/2  6 

“l  - (tl  - V'A.  7 

iNs  - (tav  - Ts,/Rs  8 

It  should  be  noted  that  ANL  and  ANS  are  truncated  to  integer  values 
in  the  computer  program.  The  smallest  of  AN^  and  ANg  is  taken  as 
the  net  change  in  passengers  between  exits  (AN) . 

5.  The  evacuation  time  for  the  longest  time  exit  is  reduced  by  AN  • 
and  the  total  amount  of  passengers  is  reduced  by  AN. 

6.  The  evacuation  time  for  the  shortest  time  exit  is  increased  by 
AN  • Rg  and  the  total  amount  of  passengers  is  increased  by  AN. 

The  same  procedure  is  again  repeated  reducing  evacuation  time  and  pas- 
sengers for  the  exit  now  having  the  greatest  evacuation  time.  The  procedure 
is  repeated  until  AN  or  AN  equals  0.  Experience  to  date  indicates  that  a 
maximum  of  five  iterations  is  required. 


PROGRAM  DESCRIPTION 


Both  wide  and  narrow  body  programs  use  the  same  names  for  subroutines. 
The  wide  body  program  has  one  more  subroutine  (UDECK)  than  the  narrow  body 
program.  The  main  program  and  subroutine  .'ATH  contain  slightly  different 
logic  between  the  two  programs.  A complete  FORTRAN  listing  of  the  wide  body 
program  is  given  in  Appendix  1.  A listing  of  the  main  program  and  subroutine 
PATH  for  the  narrow  body  program  is  given  in  Appendix  II.  A description  of 
each  subroutine  is  given  below: 
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1.  Main  program  - The  main  program  reads  in  all  the  input  data.  It 
determines  initial  passenger  exit  selection  using  two  options  de- 
scribed in  the  paragraph  entitled  INITIAL  PASSENGER  EXIT  SELECTION. 
The  main  program  calls  subroutines  GAMF,  PATH,  and  OPTIM. 

2.  Subroutine  GAMF  - This  subroutine  is  called  by  the  main  program. 

It  acts  as  a controlling  program  to  establish  the  gemma  function 
f.t  to  input  time  segment  data.  It  returns  tables  of  time  versus 
probability  to  the  main  program.  Subroutine  GAMF  calls  subroutine 
CDTR. 

3.  Subroutine  CDTR  - This  subroutine  is  called  by  subroutine  GAMF.  It 
computes  tables  of  time  versus  probability  based  on  a gamma  function 
fit  and  returns  them  to  GAMF.  Subroutine  CDTR  calls  subroutines 
NDTR  and  DLGAM. 

4.  Subroutine  NDTR  - This  subroutine  is  called  by  subroutine  CDTR.  It 
is  used  by  subroutine  CDTR  in  the  calculation  of  tables  of  prob- 
ability versus  time. 

5.  Subroutine  DLGAM  - This  subroutine  is  called  by  subroutine  CDTR.  It 
computes  the  double  precision  natural  logarithm  of  the  ganma  func- 
tion and  is  used  by  subroutine  CDTR  in  the  calculation  of  tables  of 
probability  versus  time. 

6.  Subroutine  PATH  - This  subroutine  is  called  by  the  main  program.  It 
calculates  the  evacuation  time  for  each  individual  passenger  as  de- 
scribed in  the  paragraph  entitled  CALCULATION  OF  EVACUATION  TIME. 
Subroutine  PATH  calls  subroutines  RAN DU  and  LININ. 

7.  Subroutine  RAN DU  - This  subroutine  is  called  by  subroutine  PATH. 

It  generates  a random  probability  between  0 and  1 for  each  passenger 
and  open  exit  door  which  is  returned  to  subroutine  PATH. 

8.  Subroutine  LININ  - This  subroutine  is  called  by  subroutine  PATH.  It 

uses  the  random  probability  generated  by  subroutine  RANDU  to  linearly 
interpolate  the  tables  of  probability  versus  time  generated  by  sub- 
routine GAMF.  It  returns  times  T, , T_,  T,  T. , and  T_  to  subroutine 
PATH.  1 2 3 4 5 

9.  Subroutine  OPTIM  - This  subroutine  is  called  by  the  main  program. 

It  changes  the  exit  route  of  certain  passengers  as  described  in  the 
paragraph  entitled  PASSENGER  EXIT  REDISTRIBUTION. 

10.  Subroutine  UDECX  - This  subroutine  is  only  in  the  wide  body  program. 
It  is  called  by  the  main  program  if  the  aircraft  has  an  upper  deck. 
This  subroutine  simulates  evacuation  from  an  upper  or  lower  deck 
using  the  procedure  described  in  the  paragraph  entitled  upper  deck 
wide  body  model. 
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PROGRAM  INPUT 


i 


WIDE  BODY  PROGRAM 


CARD  A FORMAT  (415) 


INPUT 

QUANTITY 

CARD  COLUMN 
NUMBERS  ^ 

i 

DESCRIPTION 

NRUN 

1-5 

No. 

of  simulations  desired 

NEXIT 

6-10 

No. 

-of  main  deck  aircraft  exits  <_ 

NOPT 

11-15 

* 0 

Passenger  evacuates  through 
nearest  open  exit 

= 1 

Passenger  evacuates  through 
assigned  exit 

(See  paragraph  entitled  initial 
passenger  exit  selection) 

NDECK 

16-20 

= 0 

No  upper  or  lower  deck 

1 Upper  or  lower  deck 


CARDS  B FORMAT  (915) 


These  cards  must  be  repeated  NEXIT  times.  The  cards  must  be  in  the  exit 
number  order  shown  in  Figure  4.  The  aircraft  seats  must  be  numbered  for  each 
exit  as  shown  in  Figure  1. 


INPUT 

QUANTITY 

CARD  COLUMN 
NUMBERS 

DESCRIPTION 

I ROWE 

1-5 

Exit  row  No.  (See  Figure  1) 

NRF 

6-10 

First  row  No.  of  passenger  seats 
assigned  to  exit  (See  Figure  1) 

NRL 

11-15 

Last  row  No.  of  passenger  seats 
assigned  to  exit  (See  Figure  1) 

NC0L1 

16-20 

Column  number  of  passenger  seats  near- 
est aisle  assigned  to  left  side  exit. 
Lowest  column  number  of  passenger  seats 
assigned  to  right  side  exit  (See  Figure 
1) 

NCOL2 

21-25 

Last  column  number  of  passenger  seats 
assigned  to  left  side  exit.  Column 
number  of  passenger  seats  nearest  aisle 
assigned  to  right  side  exit  (See 
Figure  1) 

NSIDE 

26-30 

= 0 Left  side  exit 

= 1 Right  side  exit 

NC0L3 

31-35 

Last  column  number  of  passenger  seats 
assigned  to  right  side  exit  (See 
Figure  1) 

NOPEN 

36-40 

= 0 Exit  is  open 

= 1 Exit  is  closed 

NEMP 

41-45 

No.  of  empty  passenger  seats  in  exit 

section. 


8 


CARDS  C FORMAT  (1615) 


The  I,  J values  of  empty  seats  are  input  8 pairs  to  a card.  If  there  are 
more  than  8 empty  seats  to  an  exit  continue  on  the  next  card.  A new  card 
must  be  started  for  each  exit  with  empty  seats.  Do  not  input  any  cards  for 
exits  where  NEMP  = 0. 


INPUT  CARD  COLUMN 

QUANTITY  NUMBERS  DESCRIPTION 


IV 


1-5,  11-15,  21,  I or  row  number  of  empty  seat  in 
25,  ...  section 


JV 


6-10,  16-20,  J or  column  number  of  empty  seat  in 

26-30,  ...  section 
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CARDS  D FORMAT  (5F  10.0) 


Five  sets  of  cards  are  input  for  each  open  exit  following  the  order  of 


Figure  4. 

CARD  D1 

INPUT 

CARD  COLUMN 

QUANTITY 

NUMBERS 

DESCRIPTION 

XMN1 

1-10 

Mean  value  of  time  data  for  time  Tl. 

SD1 

11-20 

Standard  deviation  of  time  data  for 
time  Tl. 

XU1 

21-30 

Maximum  value  of  time  data  for  time  Tl 

XINT1 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XUl/XINTl. 
No.  of  points  <_  50. 

XL1 

41-50 

Minimum  value  of  time  data  for  time  Tl 

CARD  D2 

XMN2 

1-10 

Mean  value  of  time  data  for  time  T2. 

SD2 

11-20 

Standard  deviation  of  time  data  for 
time  T2. 

XU  2 

21-30 

Maximum  value  of  time  data  for  time  T2 

XINT2 

31-40 

Table  interval . No.  of  points  in  time 
versus  probability  table  is  XU2/XINT2. 
No.  of  points  50. 

XL2 

41-50 

Minimum  value  of  time  data  for  time  T2 

CARD  D3 

XMN3 

1-10 

Mean  value  of  time  data  for  time  T3. 

SD3 

11-20 

Standard  deviation  of  time  data  for 
time  T3. 

XU  3 


21-30 


Maximum  value  of  time  data  for  time  T3 


INPUT 

QUANTITY 


CARO  COLUMN 
NUMBERS 


DESCRIPTION 


CARD  D3  (Continued) 


XINT3 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU3/XINT3. 
No.  of  points  <_  50. 

XL  3 

41-50 

Minimum  value  of  time  data  for  time  T3 
CARD  D4 

XMN4 

1-10 

Mean  value  of  time  data  for  time  T4. 

SD4 

11-20 

Standard  deviation  of  time  data  for 
time  T4. 

XU4 

21-30 

Maximum  value  of  time  data  for  time  T4 

XINT4 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU4/XINT4 . 
No.  of  points  < 50. 

XL4 

41-50 

Minimum  value  of  time  data  for  time  T4 
CARD  D5 

XMN5 

1-10 

Mean  value  of  time  data  for  time  T5. 

SD5 

11-20 

Standard  deviation  of  time  data  for 
time  T5. 

XU5 

21-30 

Maximum  value  of  time  data  for  time  T5 . 

XI  NTS 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU5/XINT5 . 
No.  of  points  _<  50. 

XL5 

41-50 

Minimum  value  of  time  data  for  time  T5. 

( 
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CAROS  E 


Input  only 

if  NDECK  ■ 1. 

CARD 

El 

FORMAT  (715) 

INPUT 

CARD  COLUMN 

QUANTITY 

NUMBERS 

DESCRIPTION 

I ROWE 

1-5 

Exit  row  No.  (See  Figure  3) 

NRF 

6-10 

First  row  of  passengers  in  upper  deck 
(See  Figure  3) 

NRL 

11-15 

Last  row  of  passengers  in  upper  deck 
(See  Figure  3) 

NC0L1 

16-20 

Column  No.  of  seats  on  left  of  aisle 
(See  Figure  3) 

NC0L2 

21-25 

Largest  column  No.  of  seats  in  upper 
deck  (See  Figure  3) 

NOPEN 

26-30 

« 0 Upper  deck  exit  open 

= 1 Upper  deck  exit  closed 

NEMP 

31-35 

No.  of  empty  seats  in  upper  deck 

CARD 

E2 

FORMAT  (715) 

IU 

1-5 

Row  No.  of  entrance  to  staircase  in 
upper  deck 

JU 

6-10 

Column  No.  of  entrance  to  staircase  in 
upper  deck 

IL 

11-15 

Row  No.  of  staircase  exit  on  main  deck 

JL 

16-20 

Column  No.  of  staircase  exit  on  main 
deck 

NOUT 

21-25 

Exit  No.  on  main  deck  that  upper  deck 
passengers  evacuate  through 

NNOUT 

26-30 

No.  of  passengers  on  upper  deck  that 

use  staircase 

No.  of  passengers  on  upper  deck  that  use 


NUOUT 


31-35 


CARD  E3  FORMAT  (1615) 


The  I,J  values  of  passengers  that  use  the  upper  deck  exit  are  input  eight 
pairs  to  a card.  If  there  are  more  than  eight  passengers  assigned  to  the 
upper  deck  exit  continue  on  the  next  card.  Do  not  input  if  NUOUT  = 0. 


INPUT 

CARD  COLUMN 

QUANTITY 

NUMBERS 

DESCRIPTION 

IUD 

1-5,  11-15, 

I or  row  number  of  passenger  assigned 

• 

21-25,  ... 

to  upper  deck  exit 

JUD 

6-10,  16-20, 
26-30,  ... 

J or  column  number  of  passenger 
assigned  to  upper  deck  exit 

CARD  E4  FORMAT  (1615) 


The  I,J  values  of  empty  seats  are  input  eight  pairs  to  a card.  If  there  are 
more  than  eight  empty  seats  in  the  upper  deck  continue  on  the  next  card. 

Do  not  input  if  NEMP  - 0. 


IV 

1-5,  11-15, 
21-25,  ... 

I or  row  number  of  empty  seat  in  upper 
deck 

JV 

6-10,  16-20, 
26-30,  ... 

J or  colunn  nunber  of  empty  seat  in 
upper  deck 

CARDS 

E5  FORMAT  (5F10.0) 

CARD  E5.1 

XMN1 

1-10 

Mean  value  of  time  data  for  time  T1 

SD1 

11-20 

Standard  deviation  of  time  data  for 
time  T1 

XU1 

21-30 

Maximum  value  of  time  data  for  time  T1 

XINT1 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU1/XINT1. 
No.  of  points  <_  50. 

XL1 

41-50 

Minimum  value  of  time  data  for  time  T1 

INPUT 

QUANTITY 


CARD  COLUMN 
NUMBERS 


CARD  E5.2 


\ / 


DESCRIPTION 


XMN2 

1-10 

Mean  value  of  time  data  for  time  T2  on 
upper  deck 

SD2 

11-20 

Standard  deviation  of  time  data  for 
time  T2 

XU  2 

21-30 

Maximum  value  of  time  data  for  time  T2 

XINT2 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU2/XINT2 . 
No.  of  points  <_  50. 

XL2 

41-50 

Minimum  value  of  time  data  for  time  T2 

3 not  input  cards 

E5. 3-5.5 

if  NOPEN  = 

1. 

CARD 

E5.3 

XMN3 

1-10 

Mean  -value  of  time  data  for  time  T3  of 
upper  deck  exit 

SD3 

11-20 

Standard  deviation  of  time  data  for 
time  T3 

XU  3 

21-30 

Maximum  value  of  time  data  for  time  T3 

XINT3 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU3/XINT3. 
No.  of  points  <_  50. 

XL  3 

41-50- 

Minimum  value  of  time  data  for  time  T3 

CARD 

E5.4 

XMN4 

1-10 

Mean  value  of  time  data  for  time  T4  of 
upper  deck 

• 

SD4 

11-20 

Standard  deviation  of  time  data  for 
time  T4 

• 

XU4 

21-30 

Maximum  value  of  time  data  for  time  T4 

XINT4 

31-40 

Table  interval . No.  of  points  in  time 
versus  probability  table  is  XU4/XINT4 . 
No.  of  points  <_  50. 

XL4 

41-50 

Minimum  value  of  time  data  for  time  T4 
14 

> 


CARD  E5.5 


INPUT 

CARD  COLUMN 

. 

QUANTITY 

NUMBERS 

DESCRIPTION 

XMN5 

1-10 

Mean  value  of  time  data  for  time  T5  of 
upper  deck 

S05 

11-20 

Standard  deviation  of  time  data  for 
time  T5 

• 

XU5 

21-30 

Maximum  value  of  time  data  for  time  T5 

XI  NTS 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU5/XINT5. 

No.  of  points  <_  50. 


XL5 

41-50 

Minimum  value  of  time  data  for  time  T5 

CARD  E5.6 

XMN6 

1-10 

Mean  value  of  time  data  for  time  T6  of 
staircase 

SD6 

11-20 

Standard  deviation  of  time  data  for 
time  T6 

XU6 

21-30 

Maximum  value  of  time  data  for  time  T6 

XINT6 

31-40 

Table  interval.  No.  of  points  in  time 
versus  probability  table  is  XU6/XINT6 . 
No.  of  points  <_  50. 

XL6 

41-50 

Minimum  value  of  time  data  for  time  T6 

I 


PROGRAM  INPUT 


NARROW  BODY  PROGRAM 


CARD 

A FORMAT  (415) 

INPUT 

CARD  COLUMN 

QUANTITY 

NUMBERS 

DESCRIPTION 

NHUN 

1-5 

No.  of  simulations  desired 

NEXIT 

6-10 

No.  of  aircraft  exits  <_  10 

NOPT 

11-15 

* 0 Passenger  evacuates  through 
nearest  open  exit 

= 1 Passenger  evacuates  through 
assigned  exit 

(See  paragraph  entitled  initial 
passenger  exit  selection) 

JEMP 

16-20 

Column  number  of  aisle  (See  Figure 
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CAROS  B FORMAT  (815) 


These  cards  oust  be  repeated  HEXIT  times.  The  cards  must  be  in  the  exit 
number  order  shown  in  Figure  4.  The  aircraft  seats  must  be  numbered  for  each 
exit  as  shown  in  Figure  2. 


INPUT 

CARD  COLUMN 

QUANTITY 

NUMBERS 

DESCRIPTION 

I ROWE 

1-5 

Exit  row  No.  (See  Figure  2) 

NRF 

6-10 

First  row  No.  of  passenger  seats 
assigned  to  exit  (See  Figure  2) 

NRL 

11-15 

Last  row  No.  of  passenger  seats 
assigned  to  exit  (See  Figure  2) 

NC0L1 

16-20 

Column  number  of  passenger  seats  near 
est  aisle  assigned  to  left  side  exit. 
Last  column  number  of  passenger  seats 
assigned  to  right  side  exit  (See 
Figure  2) 

NS  IDE 

21-25 

* 0 Left  side  exit 
■ 1 Right  side  exit 

NSEAT 

26-30 

■ 0 No  seats  in  exit  row 
- 1 Seats  in  exit  row 

NOPEN 

31-35 

■ 0 Exit  is  open 

■ 1 Exit  is  closed 


NEMP 


41-45 


No.  of  empty  passenger  seats  in  exit 
section 


SAMPLE  CASE  INPUT 


WIDE  BODY  PROGRAM  - Assume  an  evacuation  test  has  occurred  on  the  hypotheti- 
cal wide  body  configuration  shorn  in  Figure  5.  The  numbers  in  the  seats 
indicate  the  exit  that  passengers  used  in  the  test.  Since  we  wish  to  dupli- 
cate the  test  option  1 will  be  used  in  the  program.  For  the  purpose  of  this 
sample  case,  we  will  do  only  one  run  although  a number  such  as  ten  might  be 
better  if  we  wished  to  actually  compare  test  and  program  data.  For  card  A: 

• 

NRUN  - 1 NEXIT  » 6 NOPT  - 1 NDECK  * 0 


We  wish  now  to  assign  the  passengers  who  used  exit  1 to  that  exit.  For 
card  Bl: 


IROWE 

= 3 

NRF  * 1 

NRL  * 9 

NCOL1  = 3 

NCOL1 

= 12 

NS IDE  = 0 

NCOL  3 (blank) 

NOPEN  = 0 

NEMP  = 

■ 21 

The  only  quantity  that  needs  explanation  is  NEMP  = 21,  which  is  the 
number  of  empty  seats.  Since  we  have  assigned  the  passengers  on  the  right 
side  of  the  aircraft  to  a left  hand  exit,  the  program  logic  cannot  determine 
where  the  right  hand  aisle  is.  However,  if  we  make  column  9 empty,  we  will 
have  simulated  an  aisle.  This  totals  eight  empty  seats  since  row  3 is 
assumed  empty  by  program  logic  since  it  is  a row  exit.  For  row  8 we  wish  to 
ignore  any  passengers  who  use  an  Ci.it  different  from  1.  The  number  using 
exit  3 is  5.  Similarly  for  row  9 the  number  using  exit  3 is  8.  This  gives 
a total  of  21  empty  seats  for  exit  1. 

Although  no  passengers  are  assigned  to  exit  2,  program  logic  requires 
input  values.  We  will  input  one  empty  seat  for  exit  2.  For  card  B2: 


IROWE  = 3 

NRF  = 1 

NRL  = 1 

NCOL1  = 7 

NC0L2  = 7 

NS IDE  = 1 

NC0L3  = 7 

NOPEN  = 1 

NEMP  = 1 

We  will  later  specify  1 empty  seat  in  row  1,  column  7 for  exit  2. 
For  exit  3,  card  B3: 

IROWE  = 14  NRF  = 8 NRL  = 20  NCOL1  = 3 


NCOL2  =12  NS IDE  = 0 NC0L3  (blank)  NOPEN  = 0 

NEMP  = 28 

NEMP  = 28  is  composed  of  12  seats  for  the  right  hand  aisle,  5 seats  for 
row  8,  2 seats  for  row  9,  2 seats  for  row  19,  and  7 seats  for  row  20. 


20 


f 


For  exit  4,  card  B4: 

I ROWE  - 14  NRF  - 10  NRL  = 10  NCOL1  - 7 

NCOL2  - 7 NSIDE  - 1 NCOL3  = 7 NOPEN  = 1 

NEMP  - 1 

We  will  later  specify  1 empty  seat  in  row  10,  column  7 for  exit  4. 
For  exit  5,  card  B5: 


IROWE  - 25 
NCOL2  - 12 
NEMP  « 22 


NRF  - 19  NRL  = 30  NCOL1  =*  3 

NSIDE  = 0 NCOL3  (Blank)  NOPEN  - 0 


NEMP  = 22  is  composed  of  11  seats  for  the  right  hand  aisle,  8 seats 
for  row  19,  and  3 seats  for  row  20. 

For  exit  6,  card  B6: 

IROWE  - 25  NRF  = 20  NRL  = 20  NCOLl  = 7 

NCOL2  = 7 NSIDE  = 1 NC0L3  = 7 NOPEN  = 1 

NEMP  = 1 

We  will  later  specify  1 empty  seat  in  row  20,  column  7 for  exit  6. 

For  exit  1,  NEMP  * 21.  There  are  3 cards  Cl: 

IV  - 1,  2,  4,  5,  6,  7,  8,  9 JV  - 9,  9,  9,  9,  9,  9,  9,  9 

IV  - 8,  8,  8,  8,  8,  9,  9,  9 JV  * 1,  5,  8,  11,  12,  1,  2,  5 

IV  - 9,  9,  9,  9,  9 JV  = 6,  7,  8,  11,  12 

For  exit  2,  NEMP  * 1.  There  is  1 card  C2: 

IV  - 1 JV  = 7 

For  exit  3,  NEMP  * 28.  There  are  4 cards  C3: 

IV  - 8,  9,  10,  11,  12,  13,  15,  16  JV  = 9,  9,  9,  9,  9,  9,  9,  9 

IV  - 17,  18,  19,  20,  8,  8,  8,  8 JV  = 9,  9,  9,  9,  2,  3,  6,  7 
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IV  = 8,  9,  9,  19,  19,  20,  20,  20 


JV  = 10,  3,  10,  6,  10,  1,  2,  6 


IV  - 20,  20,  20,  20  JV  = 7,  8,  11,  12 

For  exit  4,  NEMP  = 1.  There  is  1 card  C4: 

IV  = 10  JV  = 7 


For 

exit 

5,  NEMP 

= 22 

. There 

are 

3 cards 

C5 : 

IV  = 

19, 

20, 

21, 

22, 

23, 

24, 

26, 

27 

JV  * 

9, 

9, 

9,  9,  9,  9,  9, 

9 

IV  = 

28, 

29, 

30, 

19, 

19, 

19, 

19, 

19 

JV  = 

9, 

9, 

9,  1,  2,  3,  5, 

7 

IV  - 

19, 

19, 

19, 

20, 

20, 

20 

JV  = 

8, 

11 

, 12,  3,  5,  10 

For 

exit 

6,  NEMP 

= 1. 

There 

is  1 

card 

C6 

j 

IV  * 20  JV  = 7 

Assume  that  the  following  data  have  been  obtained  for  each  exit  from 
the  evacuation  test  or  from  similar  small  scale  tests.  Note  that  input  are 


required  only  for  each 

open  exit. 

All  times 

are  in  seconds. 

EXIT 

1 

STANDARD 

MEAN 

DEVIATION 

MAXIMUM 

MINIMUM 

T1 

.6 

.6 

3. 

.3 

T2 

1. 

1. 

4. 

.5 

T3 

8.  < > 

2. 

18. 

6. 

T4 

.75 

.5 

2. 

.5 

T5 

3. 

1. 

EXIT 

3 

5. 

1. 

STANDARD 

MEAN 

DEVIATION 

MAXIMUM 

MINIMUM 

T1 

.6 

.6 

3. 

.3 

T2 

1. 

1. 

4. 

.5 

T3 

8. 

2. 

18. 

6. 
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EXIT  3 (Cont'd.) 


STANDARD 

DEVIATION 


EXIT  5 


STANDARD 

DEVIATION 


MAXIMUM 

2. 

7. 

MAXIMUM 

3. 

4. 
18. 

2. 

5. 


MINIMUM 

.5 

1.5 

MINIMUM 

.3 

.5 

6. 

.5 

1. 


Let  us  use  the  following  intervals  for  the  tine  tables  for  each  exit: 


EXIT  1 
EXIT  3 
EXIT  5 

For  cards 


.If  .1#  1. f .If  .2 
.If  .If  I. f .If  .0 
•If  .If  1. f .If  .2 


XMN1 

- .6 

SDl 

= 

.6 

XU1 

SB 

3. 

XINT1 

3 

.1 

XL1 

.3 

XMN2 

- 1. 

SD1 

s 

1. 

XU2 

SC 

4. 

XINT2 

e 

.1 

XL2 

* 

.5 

XMN3 

- 8. 

SD3 

B 

2. 

XU  3 

AS 

18. 

XINT3 

3 

1. 

XL  3 

S 

6. 

XMN4 

- .75 

SD4 

B 

.5 

XU4 

SB 

2. 

XINT4 

B 

.1 

XL4 

B 

.5 

XMN5 

- 3. 

SD5 

B 

1. 

XU5 

SB 

5. 

XINT5 

B 

.2 

XL5 

B 

1. 

For  cards 

D2: 

XMNl 

■ .6 

SDl 

B 

.6 

XU1 

3 

.3 

XINT1 

B 

.1 

XL1 

B 

.3 

XMN2 

■ .1 

SD2 

B 

1. 

XU2 

m 

4. 

XINT2 

B 

.1 

XL2 

B 

.5 

XMN3 

■ 8. 

SD3 

B 

2. 

XU  3 

SB 

18. 

XINT3 

B 

1. 

XL  3 

B 

6. 

XMN4 

- .75 

SD4 

m 

.5 

XU4 

m 

2. 

XINT4 

B 

.1 

XL4 

B 

.5 

XMN5 

- 4.5 

SD5 

m 

2. 

XU  5 

m 

7. 

XI  NTS 

B 

.5 

XL5 

S 

1.5 

For  cards  D3: 


XMN1 

■ 

.6 

SD1 

EE 

.6 

XUl 

m 

3. 

XINT1 

* 

.1 

XL1 

“ 

.3 

XMN2 

at 

1. 

SD2 

S 

1. 

XU2 

m 

4. 

XINT2 

ac 

.1 

XL  2 

= 

.5 

XMN3 

s 

8. 

SD3 

- 

2. 

XU  3 

m 

18. 

XINT3 

= 

l. 

XL  3 

m 

6. 

XMN4 

as 

.75 

SD4 

= 

.5 

XU4 

B 

2. 

XINT4 

s 

.1 

XL4 

cs 

.5 

XMN5 

3. 

SD5 

= 

1. 

XUS 

S 

5. 

XINT5 

= 

.2 

XL5 

B 

1. 

A copy  of  the  FORTRAN  coding  sheet  for  this  case  is  shown  in  Figure  6. 

UPPER  DECK  OPTION  - Assume  now  that  we  wish  to  add  an  upper  deck  to  the 
wide  body  case  previously  discussed.  Assume  that  the  upper  deck  configura- 
tion is  given  by  Figure  3.  Let  us  assign  the  eight  passengers  at  the  rear 
of  the  upper  deck  to  the  upper  deck  exit  and  the  remaining  eight  passengers 
to  the  staircase.  When  they  reach  the  main  deck  they  will  be  assigned  to 
exit  No.  1 to  evacuate  through.  Card  A must  be  changed.  Cards  E are  then 
inserted  after  cards  D for  the  main  deck. 

For  card  A: 

NRUN  « 1 NEXIT  = 6 NOPT  - 1 NDECK  * 1 

For  card  El: 

I HOWE  » 6 NRF  = 2 NRL  = 5 NC0L1  * 2 NCOL2  - 5 

NOPEN  * 0 NEMP  - 0 

Assume  the  staircase  is  at  I 8 3,  J « 6 on  the  matin  deck.  For  card 

E2: 


IU  - 1 JU-3  IL  - 3 JL  - 6 NOUT  - 1 

NNOUT  = 8 NUOUT  = 8 
For  card  E3: 


IUD  - 4,  4,  4,  4,  5,  5,  5,  5 
JUD  - 1,  2,  4,  5,  1,  2,  4,  5 
Since  NEMP  - 0,  there  is  no  card  E4. 

Assume  the  following  data  has  been  obtained  from  evacuation  tests : 
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UPPER  DECK  EXIT 


MEAN 

STANDARD 

DEVIATION  MAXIMUM 

MINIMUM 

T1  .6 

.6 

3. 

.3 

T2  1. 

1. 

4. 

.5 

T3  12. 

4. 

20. 

6. 

T4  1.5 

.5 

4. 

.75 

T5  5. 

1. 

10. 

2. 

STAIRCASE 

T6  4. 

1. 

8. 

1. 

Let  us  use  the 

following  intervals  for  the  time 

tables t .1, 

. 1#  .2/  .2 , .5,  .2 

For  cards 

E5: 

XMNl  - .6 

SD1 

- .6 

XU1  - 3. 

XINT1  - .1 

XL1  - .3 

XMN2  - 1. 

SD2 

- 1. 

XU2  - 4. 

XINT2  - .1 

XL2  - .5 

XMN3  - 12. 

SD3 

* 4. 

XU 3 - 20. 

XINT3  - 2. 

XL 3 - 6. 

XMN4  - 1.5 

SD4 

= .5 

XU4  * 4. 

XINT4  - .2 

XL4  = .75 

XMN5  - 5. 

SD5 

“ 1. 

XU5  - 10. 

XINT5  - .5 

XL5  = 2. 

XMN6  - 4. 

SD6 

- 1. 

XU6  - 8. 

XINT6  - .2 

XL6  = 1. 

The  FORTRAN  coding  for  this  case  is  shown  in  Figure  6 after  the  coding 
for  the  main  deck  only  case. 

NARROW  BODY  PROGRAM  - Assume  we  wish  to  simulate  an  emergency  evacuation  in 
the  narrow  body  aircraft  shown  in  Figure  7.  In  this  case  we  will  let  pas- 
sengers evacuate  through  the  nearest  open  exit.  Seating  assignment  to  an 
exit  is  not  important  except  that  all  the  seats  in  the  aircraft  must  be 
assigned  to  exits. 

Let  us  simulate 

ten  evacuations. 

For  card  At 

NRUN 

- 10 

NEXIT  - 8 NOPT  - 0 

JEMP  = 4 

nn 


! 


Let  us  assign  the  first  class  passengers  to  exits  1 and  2.  For  card 

Bit 

IROWE  - 1 NRF  - 2 NRL  - 8 NCOL1  - 3 


NS IDE  - 0 

NSEAT  - 0 

NOPEN  - 0 

NEMP  - 7 

NEMP  « 7 will  be  used  to 
first  class  seating. 

delete  passengers  in  column 

1 thus  simulating  the 

For  card  B2: 

IROWE  - 1 

NRF  * 2 

NRL  * 8 

NCOL1  - 6 

NS IDE  - 1 

NSEAT  - 0 

NOPEN  - 1 

NEMP  - 0 

It  is  not  necessary  to  delete  any  passengers  on  this  side  of  the  aircraft 
since  program  logic  will  take  columns  5 and  6 only  because  NC0L1  - 6.  For 
the  left  side  exit  NCOL1  ■ 3 and  the  program  logic  will  take  columns  1,  2, 
and  3.  If  NCOL1  * 2 had  been  input  there  would  have  been  a mismatch  with 
the  aisle  column  which  is  input  as  JEMP  - 4. 

For  card  B3: 

IROWE  - 12 

NRF  - 9 

NRL  - 14 

NCOL1  - 3 

NS IDE  - 0 

NSEAT  = 1 

NOPEN  - 1 

NEMP  - 0 

Notice  there  are  seats  in  the  exit  row  for  this  exit 

SO  NSEAT  - 1. 

For  card  B4 : 

IROWE  - 12 

NRF  « 9 

NRL  - 14 

NCOL1  - 7 

NS IDE  - 1 

NSEAT  - 1 

NOPEN  - 0 

NEMP  - 0 

For  card  B5: 

IROWE  - 15 

NRF  - 15 

NRL  - 20 

NCOL1  - 3 

NS IDE  - 0 

NSEAT  - 1 

NOPEN  - 1 

NEMP  - 0 

For  card  B6 : 

IROWE  - 15 

NRF  - 15 

NFL  - 20 

NCOL1  - 7 

NS IDE  - 1 

NSEAT  - 1 

NOPEN  - 0 

NEMP  - 0 
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For  card  B7: 


IROWE  - 27  NRF  - 21 

NS IDE  - 0 NSEAT  - 0 

For  card  B8: 

IROWE  - 27  NRF  - 21 

NS IDE  - 1 NSEAT  - 0 

For  card  Cl: 

IV  - 2,  3,  4,  5,  6,  7,  6 
There  are  no  more  C cards. 


NRL  - 26 
NOPEN  - 0 

NRL  - 26 
NOPEN  - 1 


NCOLl  - 3 
NEMP  - 0 

NCOLl  - 7 
NEMP  - 0 


OV  ■ 1,  1,  1,  1,  1,  1,  1 


Assume  that  the  following  data  are  available  from  previous  testing. 
There  are  four  sets  of  this  data,  one  set  for  each  open  exit.  All  times 
are  in  seconds. 


MEAN 

EXIT  1 

STANDARD 

DEVIATION 

MAXIMUM 

MINIMUM 

T1 

.6 

.6 

3. 

.3 

T2 

1. 

1. 

4. 

.5 

T3 

10. 

2. 

20. 

6. 

T4 

1. 

1. 

3. 

.5 

T5 

2.5 

1. 

4. 

1. 

• 

EXIT  4 

* 

STANDARD 

MEAN 

DEVIATION 

MAXIMUM 

MINIMUM 

T1 

.6 

.6 

3. 

.3 

T2 

1. 

1. 

4. 

.5 

T3 

6. 

2. 

12. 

6. 

T4 

1.5 

1.5 

4. 

.5 
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EXIT  4 Oont'd 


MEAN 

STANDARD 

DEVIATION 

MAXIMUM 

MINIMUM 

T5 

4. 

1. 

6. 

3. 

EXIT  6 


T1 

MEAN 

.6 

STANDARD 

DEVIATION 

.6 

MAXIMUM 

3. 

MINIMUM 

.3 

T2 

1. 

1. 

4. 

.5 

T3 

6. 

2. 

12. 

6. 

T4 

1.5 

1.5 

4. 

.5 

T5 

4. 

1. 

6. 

3. 

T1 

MEAN 

.6 

STANDARD 

DEVIATION 

.6 

EXIT  7 

MAXIMUM 

3. 

MINIMUM 

.3 

T2 

1. 

1. 

4. 

.5 

T3 

10. 

2. 

20. 

6. 

T4 

1. 

1. 

3. 

.5 

T5 

2.5 

1. 

4. 

1. 

Let  us 

use  the 

following  intervals 

for  the 

time  tables 

for  each  exit: 

EXIT  1 .1,  . 1,  2.,  .1,  .2 
EXIT  3 .1,  .1,  1.,  .1,  .2 
EXIT  6 .1,  .1,  1.,  .1,  .2 
EXIT  7 .1,  .1,  2.,  .1,  .2 


For  cards 


XMN1 

- .6 

SD1 

m 

.6 

XUl 

as 

3. 

XXNT1 

m 

.1 

XL1 

XMN2 

- 1. 

S02 

m 

l. 

XU2 

m 

4. 

XINT2 

m 

.1 

XL2 

XMN3 

- 10. 

SD3 

m 

2. 

XU  3 

m 

20. 

XINT3 

m 

2. 

XL  3 

XMN4 

- 1. 

SD4 

m 

1. 

XU4 

m 

3. 

XINT4 

m 

.1 

XU 

XMN5 

- 2.5 

S05 

m 

1. 

XU5 

m 

4. 

XINT5 

m 

.2 

XL5 

For  cards 

D2: 

XMN1 

- .6 

SOI 

m 

.6 

XUl 

m 

3. 

XINT1 

m 

.1 

XLl 

XMN2 

- 1. 

SD2 

m 

1. 

XU  2 

m 

4. 

XINT2 

m 

.1 

XL2 

XMN3 

- 6. 

S03 

m 

2. 

XU3 

m 

12. 

XI  NT  3 

m 

1. 

XL  3 

XMN4 

- 1.5 

SD4 

m 

1.5 

XU4 

m 

4. 

XINT4 

m 

.1 

XU 

XMN5 

- 4. 

SD5 

m 

1. 

XUS 

m 

6. 

XI  NTS 

m 

.2 

XL5 

For  cards 

03 1 

XMN1 

« .6 

SDl 

m 

.6 

XUl 

m 

3. 

XINT2 

m 

.1 

XLl 

XMI2 

- 1. 

SD2 

m 

1. 

XU2 

m 

4. 

XINT2 

m 

.1 

XL2 

XMN3 

- 6. 

SD3 

m 

2. 

XU  3 

m 

12. 

XINT3 

m 

1. 

XL  3 

XMN4 

» 1.5 

S04 

m 

1.5 

XU4 

m 

4. 

XXNT4 

m 

.1 

XU 

For  cards 

04: 

XMN1 

- .6 

SDl 

m 

.6 

XUl 

m 

3. 

XINTl 

m 

.1 

XLl 

XMN2 

- 1. 

SD2 

m 

1. 

XU2 

m 

4. 

XINT2 

m 

.1 

XL2 

XMN3 

- 10. 

SD3 

m 

2. 

XU  3 

m 

20. 

XINT3 

m 

2. 

XL  3 

XMN4 

- 1. 

SD4 

m 

1. 

XU4 

m 

3. 

XINT4 

m 

.1 

XU 

XMN5 

- 2.5 

SOS 

m 

1. 

XUS 

m 

4. 

XINT5 

m 

.2 

XL5 

A copy  of  the  FORTRAN  coding  sheet  for  this  case  is  shown  in  Figure  8. 


PROGRAM  OUTPUT 


A copy  of  the  program  output  for  the  wide  body  sample  case  with  upper 
deck  is  shown  in  Figure  9.  The  first  page  of  output  consists  of  a print  out 
of  the  input  aircraft  configuration  data  in  the  same  format  as  the  input. 

This  allows  the  user  an  easy  check  to  make  sure  input  data  is  correct.  The 
aircraft  configuration  data  are  followed  by  the  time  segment  input  data. 

The  exit  appropriate  to  the  time  data  is  identified  in  the  last  column  of 
this  output.  Next  are  the  input  data  for  the  upper  deck  if  this  option  is 
used  in  the  program. 

The  first  output  consists  of  a list  of  passengers  identifying  seat  row, 
seat  column,  exit  assigned,  and  exit  out.  For  each  passenger  his  time  to 
the  exit  door,  door  delay  time,  and  evacuation  time  are  given.  By  comparing 
time  to  the  exit  door  and  door  delay  time  the  user  can  determine  if  a bottle- 
neck is  occurring  at  the  door.  If  time  to  the  exit  door  is  greater  than  door 
delay  time  for  any  passenger,  this  indicates  that  the  door  is  not  always  a 
bottleneck.  A new  page  is  started  for  each  exit  that  is  used  in  the  evacua- 
tion. The  first  exit  listed  is  the  upper  deck  exit  if  the  upper  deck  option 
is  used.  This  exit  number  is  set  equal  to  NEXIT  + 1.  This,  passengers  on 
the  upper  deck  that  use  a main  deck  exit  may  be  identified  by  looking  for 
NEXIT  + 1 in  their  exit  assigned  column.  The  remaining  exit  data  is  in  the 
order  of  main  deck  exit  numbers. 

The  next  page  of  output  starts  out  with  a summary  of  the  total  number  of 
passengers  that  used  an  exit,  the  total  evacuation  time  for  that  exit,  and 
the  exit's  number  for  the  initial  evacuation.  Passengers  are  then  redistrib- 
uted to  exits  with  shorter  lines  and  the  next  line  of  output  indicates  where 
passengers  were  taken  from  and  where  they  went.  The  new  passenger  totals  and 
evacuation  times  for  all  the  exits  are  then  listed.  This  is  repeated  until 
logic  is  not  able  to  redistribute  any  more  passengers  between  exits. 

The  narrow  body  program  has  the  same  output  as  the  wide  body  program 
except  there  is  no  upper  deck  output  for  the  narrow  body  program.  A copy  of 
the  program  output  for  one  of  the  ten  simulations  of  the  narrow  body  sample 
case  is  shown  in  Figure  10.  All  times  are  in  seconds  in  these  programs.  The 
CPU  time  for  five  evacuations  is  approximately  one  minute  on  the  IMB  370/155 
for  this  case.  CPU  time  has  ranged  up  to  90  minutes  for  100  evacuations  of  a 
389  passenger  wide  body  aircraft. 
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SUMMARY 


Mathematical  models  for  wide  and  narrow  body  aircraft  emergency  evacua- 
tions have  been  developed.  Fortran  computer  programs  have  been  developed 
from  the  mathematical  models.  Several  cases  have  been  run  using  all  program 
options  producing  reasonable  results  for  the  wide  and  narrow  body  programs. 
The  programs  require  evacuation  path  time  segment  input  data.  Some  small 
scale  testing  will  be  required  to  generate  a valid  range  of  appropriate  time 
segment  input  data.  Correlation  between  full  scale  evacuation  tests  and 
program  output  should  then  be  performed  to  provide  program  validation. 


I = 
NRF 

NRL 

IROWE 


3 

1 

2 

3 

4 

5 

6 


NCOL1  NCOL2 
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Fig.  3 - Upper  deck  model 
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